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Goal 


• Goal: Design an optimal powered descent trajectory on-board the 
spacecraft in order to softly land on an asteroid. 

■ Algorithm would need to be autonomous, reliable, robust and efficient. 

■ Designing on-board facilitates easy change of parameters. 

• Convex optimization is efficient and reliable. 

■ Guarantees global minimum in finite number of steps, if problem is feasible. 

■ Especially true in subclasses, such as Second Order Cone Programming. 

• Can convex optimization be used to design the asteroid powered 
descent trajectory? 


Original Problem Statement 


Asteroid powered descent fuel optimal problem is nonlinear and 


nonconvex. 


S.t. r = V, V = - — 
1 m 


min —m(tf) 

2uoxv — uoxr — uox (w x f) + VU (r), m = — - 


T < 

mm _ 


T 


T 


< T 

— -‘-max-) 


r 


cosd — r T n < 0, to > mdry 


r(0) = f 0 , v (0) = vq, to (0) = m wet , f{t f ) = r f , v(t f ) = v f , t f given 


Highlighted terms are not permissible for a convex optimization problem. 


Relaxation 


• Original Problem 

min —m(tf) 
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Slack Variable 


New Constraint 
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Proved the optimal solution of the relaxed problem is the optimal solution of the original. 





Successive Solution Method 


• Solve a series of convex optimization problems. 

• Equations of motion can be arranged as: 
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• A and c are evaluated using the previous solution (k-1). 

• In the (k) th iteration, dynamics are linear and time varying. 

• Iterations continue until two successive trajectories are within a set 
tolerance. 


Scaling and Discretization 


• Position, velocity, acceleration and time are scaled to keep the 
magnitudes near 1. 

■ Position scaled by smallest semi-major axis 

• Discretize the continuous equations to form equality constraints. 

• Now, the problem is in the form of a convex optimization problem. 

■ Linear equality constraints 

■ Convex inequality constraints 

• Used CVX, publically available Matlab based convex solver. 

■ A 600 second trajectory, 300 nodes ~ 3300 variables per iteration, 4 iterations 
takes 3.5 minutes to complete. 


Asteroid Assumptions 


• Two Trajectories 

■ North Pole Trajectory (NP) lands on the z axis 

■ Equatorial Trajectory (EQ) lands on the x axis 


• Three Asteroid Sizes 

■ Al: 1000 x 500 x 250 m 

■ A2: 750 x500 x250 m 

■ A3: 500 x 500 x 250 m 

• Four Periods 

■ 12 hr, 8 hr, 4 hr, 2 hr 

■ Constant spin rate along z axis. 
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Propellant Usage Parameter Sweeps 
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EQ trajectories 1000 x 500 x 250 m 
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EQ Optimal Trajectory 


• Minimum Propellant flight time 512 sec. 

• 1000 x 500 x 250 m 
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Conclusions 


• Powered descent trajectory can be formulated as a convex 
optimization problem. 

• Successive solution method is the key to handling nonlinear gravity 
model. 

• Handles a wide range of parameters successfully. 

• Viable algorithm for rapidly designing asteroid powered descent 
trajectories on-board the spacecraft. 


Back-up 


Coordinates 



Use an asteroid centered fixed coordinate system to determine the spacecraft location, 
gravity and equations of motion. 

• x aligned with the largest axis, z aligned with the smallest axis, y completes the orthogonality. 

• Radius (r) is the distance from asteroid center to spacecraft. 

• A spacecraft longitude 

• 6 spacecraft latitude 


Gravitational Acceleration 


Equations 
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Where 

r = (rl+r% + rf) 
6 = arcsin (^f) 
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Spacecraft Parameters 


• Thrust 80 N - 20 N 

■ Isp 225 sec 

• Mass 1400 kg 

• Initial Conditions 
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